Astronomy & Astrophysics manuscript no. M31 


© ESO 2009 


July 3, 2009 





Dust properties in M31.I. 
Basic properties and a discussion on age-dependent dust heating 

M. Montalto 1 - 2 , S. Seitz 1 ' 2 , A. Riffeser 1 , U. Hopp 1 - 2 , C.-H. Lee 1 , R. Schonrich 1 

cs ■ 

University Observatory Munich, Scheinerstrasse 1, 81679 Munich, Germany e-mail: montalto@usm.uni-muenchen.de 
^ ■ 2 Max-Planck Institute for Extraterrestrial Physics, Giessenbachstrasse, 85748 Garching, Germany 



July 3, 2009 

ABSTRACT 



Context. 

Aims. Recent observations derived from the Spitzer Space Telescope and improvements in theoretical modeling of dust emission 
properties are used to discuss the distribution of dust and its characteristics in the closest neighbor spiral galaxy M3 1 . Together with 
T ^ ' GALEX FUV, NUV, and SDSS images we studied the age dependence of the dust heating process. 

Methods. 

Spitzer IRAC/MIPS maps of M3 1 were matched together and compared to dust emission models allowing to constrain the dust mass, 
[■ the intensity of the mean radiation field, the abundance of Polycyclic Aromatic Hydrocarbons (PAH) particles. The total infrared 

emission (TIR) was analyzed in function of UV and Optical colors and compared to predictions of models which consider the age- 
dependent dust heating. 

Results. We demonstrate that cold-dust component emission dominates the infrared spectral energy distribution of M3 1 . The mean 
intensity of the radiation field heating the dust is low (typically U < 2, where U = 1 is the value in the solar surrounding). Due to 
the lack of submillimetric measurements the dust mass (M t | llst ) is only weakly constrained by the infrared spectrum, but we derived a 
lower limit of M(j ust > 1 . 1 x 10 7 M o with a best fit value of M t | llst = 7.6 x 10 7 M Q , in good agreement with expectations from CO and HI 
measurements. We show that across the spiral-ring structure of M3 1 a fraction > 3% of the total dust mass is in PAHs. UV and optical 
colors are correlated to total infrared to far ultraviolet (TIR/FUV) ratios in ~ 670pc-sized regions overall the disk of M31, although 
' deviating from the relationship between infrared excess and ultraviolet spectral slope (referred as IRX -f} relationship) for starburst 

galaxies. In particular, redder regions show lower values of the TIR/FUV ratio for a fixed color. Considering the predictions of models 
that account for the dust-heating age dependence we derived that in 83% of the regions analyzed across the lOkpc ring more than 50% 
of the energy absorbed by the dust is rediated at A > 4000A and that dust in M31 appears mainly heated by populations a few Gyr old 
even across the star-forming ring. We also found that the attenuation is varying radially peaking near lOkpc and decreasing faster in 
the inner regions of M3 1 than in the outer regions in agreement with previous studies. We finally derived the attenuation map of M3 1 
at 6 /px resolution (~ 100 pc/px along the plane of M31). 
Conclusions. 
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H I 1- Introduction al. 2007), SDSS optical images (York et al. 2000), and Spitzer 

Space Telescope observations (Barmby et al. 2006, Gordon et 

The presence of the dust affects astrophysical observations in al 2006) in the infrared and far infrared are now avai i able for 

different ways. While in the UV-optical spectral region the dust M3 { Toge ther with the recent CO and HI observations (Nieten 

absorbs the stellar radiation field, it re-emits the absorbed energy et al 2 006, Braun et al. 2009) this new wealth of data allows a 

in the IR-FIR and sub-millimetric spectral regions. It is thus es- much more detailed exploration of the dust properties out than 

sential to determine the properties of the dust in a given system in prev i ous stud ies. Especially Spitzer data allow to study signif- 

as accurately as possible. While in our own galaxy we can reach icantly smaller sca les than in the IRAS era. 

very high spatial resolutions and accuracies, our particular loca- Moreove r, in the last years important progress has been made 

tion inside the galaxy limits our analysis. on the the oretical understanding of dust properties facilitating 

Being the closest large spiral like the Milky Way, M31 of- me derivation of new parameters from the data (Weingartner & 

fers a fascinating view of the dust distribution across the disk of Draine 2001; Draine & Li 2007; Draine et al. 2007). This work 

the galaxy and an unique opportunity to study the properties of a i s thus devoted to an analysis of recent data available for M3 1 in 

major disc galaxy in detail. Previous studies already focused on light of the new theoretical instruments. We focus our attention 

the dust in this galaxy, but new observational results and theoret- on the mean radiation field heating the dust, on the dust mass 

ical improvements require some issues to be revisited and criti- determination and on the amount and distribution of Polycyclic 

cally discussed and allow to shed new light onto still unresolved Aromatic Hydocarbon particles (PAHs) in M3 1 . 
problems. In particular GALEX UV observations (Gil de Paz et M3 j offers a nice example for studying a problem that has 

drawn a lot of attention in recent studies: what is the main source 

Send offprint requests to: heating the dust in a galaxy? Certainly the dust tends to prefer- 
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antially absorb UV photons, so naively one would of course ex- 
pect UV radiation especially from the young populations seen 
in the GALEX images of M31 to be the main source power- 
ing the observed dust emission. Indeed this idea clearly holds 
for starburst and very actively star forming galaxies (see e.g. 
Buat 1992, Meurer et al. 1995) and triggered the discovery of 
the IRX-/3 relation (Meurer et al. 1999), but in galaxies host- 
ing older stellar populations the radiation field is clearly domi- 
nated by longer wavelengths. Recent studies have demonstrated 
that in systems with low relative star formation rate old stars can 
contribute significantly to the dust heating (Cortese et al. 2008; 
Kong et al. 2004, Buat et al. 2005; Gordon et al. 2000) and 
we will see in this paper that the latter applies for M31. This 
difference is crucial if one attempts to derive accurate dust at- 
tenutations and consequently star formation rates in any galactic 
system. 

The structure of this paper is as follows: In Sect. 2 we present 
the examined data and the calibration procedure, in Sect. 3 we 
describe the dust models and their application to determine the 
dust mass, the mean intensity of the radiation field responsi- 
ble for dust heating and of the PAHs abundances in M31. In 
Sec. 4 we discuss the age dependency of dust heating, deriv- 
ing the A FU v attenuation across the spiral-ring structure of M3 1 . 
Finally in Sects. 5 and 6 we discuss and summarize our results. 

2. The data 

2.1. IR 

In this work we used the IR images of M3 1 described in Barmby 
et al. (2006) and Gordon et al. (2006), and obtained with the 
IRAC/MIPS instruments on board of the Spitzer Space Telescope 
(SST). In particular, we retrieved from the SST archive the 
3.6 pm, Spm IRAC, and the 24 pm, 70 pm and 160yum MIPS 
Basic Calibrated Data (BCDs) which were processed by the 
SST's archive pipeline (v. 14). The BCDs were then stacked to- 
gether using the MOPEX software (version 18.1). 

For IRAC observations we applied overlap correction before 
mosaicing in order to match the background between adjacent 
images and discarded the first frame of each observing run be- 
cause it is typically acquired with a lower bias level than subse- 
quent images (first frame effect). The background was evaluated 
in the external regions and subtracted from the images with the 
IRAF imsurfit task using a low order polynomial (first or second 
order depending on the map). We then applied the IRAC pho- 
tometric corrections for infinite sources as reported in the IRAC 
documentation multiplying the 3.6 pm image by 0.91 and the 
8 pm image by 0.74 in order to account for scattering of incident 
light in the array focal planes. We applied the correction for ex- 
tended sources because in this work we focus on the diffuse dust 
emission, not on point-like sources. 

For MIPS 70 pm and 160/im data, we used our customized 
software for background subtraction via fitting and subtracting 
a first order polynomial along the scan direction excluding the 
M3 1 region. The final MIPS mosaic images have a resolution of 
6" /px, 18" /px and 40" / px at 24 pm, 70 pm and 1 60 pm, respec- 
tively. 

Because our purpose was to directly compare the fluxes in 
each band on a pixel by pixel basis we had to resample the im- 
ages to the same astrometric reference system. In fact we created 
two sets of images: in the first one the IRAC 3.6 pm and 8 pm 
images were put to the reference system of the 24 pm mosaic 
(6"/px), in the second one the IRAC 3.6 pm and 8 pm and the 
MIPS 24 pm and 70 pm were put on the reference system of the 



Table 1. Percentual variation of the total integrated flux of M31 as mea- 
sured before and after the resampling process of the images described 
in the text. The third column presents the percentual variation for the 
images resampled to the 24 pm resolution, the third column for the im- 
ages resampled to 160/imi. The total integrated flux was measured in 
the same area in all the images, correspondent to an ellipse centered on 
M31, with a semi-major axis of 84', a semi-minor axis of 18', and a 
position angle of PA = 38°. 



Band (pm) 


resampled to 24 pm (%) 


resampled to 160 pm (%) 


3.6 


2 


2 


7.9 


3 


2 


24 




0.8 


70 




0.6 


160 







160 pm mosaic (40" /px). These two sets will be considered sep- 
arately in the following. To make the different instrument images 
comparable, we first matched the PSF of the input instrument to 
that one of the target instrument (MIPS 24 pm or MIPS 160 pm). 
This was done with the help of the transformation kernels pro- 
vided by K. Gordon 1 (Gordon et al. 2008). These kernels give 
the transformation between the PSFs of each couple of IRAC- 
IRAC, IRAC-MIPS and MIPS-MIPS instruments dependent on 
the dust temperature. After the convolution we resampled the im- 
ages to the same astrometric grid, in such a way that each pixel 
in a given band corresponded to the same physical region of the 
same pixel in the other bands. We estimated the uncertainty of 
the resampling step by means of a comparison of the integrated 
photometry on the whole galaxy before and after this procedure 
and found that in the worst case the relative variation in flux was 
3%, as shown in Tab. 1. The uncertainty in the background sub- 
traction was estimated evaluating the scatter in the background 
subtracted mosaics in the outermost regions and found to be neg- 
ligible. 

While creating the mosaics with MOPEX 2 (MOsaicker and 
Point source Extractor), we derived also the associated error 
map, which provides the error related to the whole reduction 
procedure performed by the pipeline. To account for the other 
sources of error (resampling, background subtraction, BCDs cal- 
ibration) we simply added a constant term equal to 10%, with 
except of the 24 pm map for which we would have otherwise ob- 
tained unrealistically large values of the error. In reality our ap- 
proach can be considered rather conservative, because as stated 
in the SST documentation the standard deviation maps account 
for the error of the pipeline in a conservative way. This is also 
demonstrated by a comparison of our own measurements of the 
integrated fluxes of M31 in each band with those derived by 
Barmby et al. (2006) and Gordon et al. (2006) as shown in Tab. 2. 
In particular, with respect to the errors quoted by Barmby et 
al. (2006) and Gordon et al. (2006), our integrated flux mea- 
surements are 0.6<x smaller, O.lcr smaller, lcr larger, 0.8<x larger 
and 0.4cr smaller than their own measurements at 3.6 pm, 8 pm, 
24pm 10pm and 160/im respectively, and our estimated error 
is 9% smaller, 30%, 54%, 36% and 50% larger than their given 
uncertainties. We conclude that our measurements are consistent 
within lcr and that our errors are in general equivalent or larger 
with respect to the values derived by Barmby et al. (2006) and 
Gordon et al. (2006). 



1 http://dirty.as.arizona.edu/ kgordon/mips/conv_psfs/conv_psfp.html 

2 http://ssc.spitzer.caltech.edu/postbcd/mopex.html 
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We finally subtracted from the 8 pm and the 24 pm maps the 
stellar continuum which was assumed to be described by the 
3.6 pm map following Helou et al. (2004): 



F° s (8 pm) = F v (8 pm) - 0.232F V (3.6 pm) 
F™(24pm) = F v (24pm) - 0.032^(3.6^) 

2.2. Ultraviolet 



(1) 
(2) 



The UV maps of M31 were retrieved directly from the Nasa 
Extragalactic Database (NED) and were obtained as part of the 
survey of nearby galaxies performed with the GALaxy Evolution 
eXplorer (GALEX) satellite (Gil de Paz et al. 2007). We re- 
trieved and used both the far ultraviolet (FUV) and the near ultra- 
violet (NUV) images. These maps have a pixel scale of 1 .5" /px, 
a total integration time of 800 sec and are given in units of counts 
per pixel per second (CPS). The background in these images can 
be assumed constant and was accurately measured by Gil de Paz 
et al. (2007) who give a mean value of (9.29 + 5.24) x 10~ 4 
CPS and 4.712 + 1.529) x 10~ 3 CPS for the FUV and NUV 
respectively. We subtracted the background from the maps and 
then converted the CPS units in Jy multiplying by the factors 
/o.fuv = 108yuJy/CPS and /o.nuv = 36/iJy/CPS, as given by 
the GALEX documentation. We then accounted for the fore- 
ground galactic extinction assuming E(B - V) = 0.062 in the 
direction of M31 as reported by Schlegel et. al. (1998) with 
A FUV = 7.9 E(B - V) and A NUV = 8.0 E(B - V) as given by 
Gil de Paz et al. (2007). Two series of maps were finally de- 
rived in order to compare UV observations with MIPS 24 pm and 
MIPS 160 pm data. Both the FUV and the NUV maps were con- 
volved to the correspondent MIPS PSFs and regridded at 6"/px 
and 40"/px. The global uncertainty of the final maps was esti- 
mated to be ~ 10% in flux for both maps. 

2.3. Optical 

We retrieved and stacked together the SDSS images of M31 ob- 
tained in the SDSS survey. In this work we made use only of 
the SDSS i band (hereafter /sdss) images of M31. First we sub- 
tracted the SOFTBIAS in every single frame. SOFTBIAS was 
originally added automatically by the SDSS pipeline to avoid 
negative pixel values. We aligned each single frame to a large 
reference image (1.5 x 4.2 square degree) using the IRAF task 
wregister. Then we reconstructed each stripe. The final mo- 
saic image contains 7 stripes (3366-2, 3367-3, 3366-3, 3367-4, 
3366-4, 3367-5, and 3366-5, from left to right in the mosaic im- 
age). We first tested for photometric variation, but resulted in no 
zero-point difference between the stripes. Nevertheless we found 
sky variations in y-direction in different stripes when combining 
them into the mosaic image. Therefore, we applied an alignment 
in y-direction for each stripe. Using our customized written soft- 
ware (mupipe/skycalc) we calculated the gradient a x y and the 
offset b in flux of the overlapping areas between two adjacent 
stripes, e.g. F3366-2 = ^3367-3 + «3367-3 xy + b 3 367-3 and applied 
the y-dependent sky correction to the whole stripe. The map was 
convolved with the Spitzer 160 pm PSF and rescaled at 40" /px 
resolution. 

Then we calibrated the images using cataloged surrounding 
stars with SDSS photometry and masked out all the sources with 
; < 17 mag as most of the point-like sources in this magnitude 
range are foreground contaminants and excluding the brightest 
M31 objects (globular clusters, red supergiants, etc.) does not 
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Fig. 1. Two color diagram used to isolate the regions with dust diffuse 
emission. Labels and dashed lines indicate different sources of emis- 
sion in function of the colors. (S/N) > 1 in all maps and each region 
corresponds to 6"/px resolution. 



affect our analysis which is focused on regions dominated by the 
dust diffuse emission. 

In order to select a un-contaminated sample of regions from 
the infrared and UV maps we used the two color diagram shown 
in Fig. 1. We corrected the 8 pm and 24 pm maps for the stellar 
continuum using Eq. 1 and Eq. 2 and worked at 6" /px resolution. 
We limited the analysis to regions with (S /N) > 1 by this exclud- 
ing from the analysis the inner bulge region (< 3' from the center 
of M31), partially the interarm regions between the inner spiral- 
arm and the outermost regions of the galaxy. While the bulk of 
the points have colors with 1 < v F" s (8pm)/v Ff(24pm) < 10 
and 0.3 < FUV/NUV < 2, and appear after visual inspection 
associated with the diffuse dust emission, some deviant points 
are visible in Fig. 1 . In our maps bright point-like sources in the 
24 pm map have v Ff(%pm) / V Ff(24pm) < 1, thus lower than 
the ratios found in the diffuse component of the dust. Different 
studies have demonstrated that the emission from point-like 
sources at 24 pm is associated with the emission from HII re- 
gions in the UV and Optical spectral regions (Calzetti et al. 2007; 
Prescott et al. 2007). Moreover, the PAHs emission at 7.7 pm is 
expected to be strongly reduced with respect to the diffuse dust 
emission in such environments (e.g. Calzetti et al. 2005; Thilker 
et al. 2007; Bendo et al. 2008). On the other hand regions with 
vF™(8pm)/vF™(24pm) > 10 or FUV/NUV < 0.3 appear to be 
associated with bright foreground stars. We thus adopt the selec- 
tion criteria for regions: 1 < v F" s (8pm)/v F™(24pm) < 10 and 
FUV/NUV > 0.3. Once using also the SDSS optical map we 
further added the selection for excluding bright foreground stars 
as described above. 

In Fig. 2 for each of the analyzed regions of M31 is shown 
the minimum signal to noise ratio (S/N) when considering all 
the above mentioned observational maps together. The (S /N) is 
in general larger than 3 along the lOkpc ring and partially in the 
the inner spiral structure. In the interarm and external regions 
the (S/N) drops below 1.5. In the bulge region the S/N is low 
because of the small stellar corrected flux at 8 pm. 
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U.b 1 l.b Z l.b b b.b 4 4.ij 

Fig. 2. Upper panel: Minimum S/N in each one of the analyzed regions of M31 once considering the FUV, NUV, isDss. 3.6yum, 8/jm, 24pm, 70pm 
and 160/jm observational maps described in Sec. 2 all together. The resolution is 6 /px. We show also the nine concentric 2kpc-wide regions used 
in Sec. 4. Lower panel: the same as the upper panel but for the 40 /px resolution maps. 



3. The theoretical models 

In this work we used and tested the theoretical models of in- 
frared emission of dust grains heated by starlight developed by 
Weingartner & Draine (2001) and Li & Draine (2001), recently 
updated by Draine & Li (2007). In the following we summarize 
their main characteristics and results. Dust grains are treated as 
a mixture of amorphous silicate and carbonaceous grains, where 
the smallest carbonaceous grains have the physical properties 
of the PAHs. Temperature distribution functions for all particles 
are calculated and emission spectra provided for dust heated by 
a single radiation field intensity. The spectral energy distribution 
for this is given by Mathis, Mezger, & Panagia (1983) (MMP83) 
scaled by the dimensionless factor U, where U — 1 corresponds 
to the mean intensity of the radiation field in the local ISM of the 
solar neighborhood. We tested the seven Milky Way dust mod- 
els proposed by Weingartner & Draine (2001), which have PAH 
abundances with 0.47%, 1.12%, 1.77%, 2.50%, 3.19%, 3.90% 
and 4.58% of the total dust mass. These models assume a dust- 
to-gas ratio Mdust/^H — 0.01. For further details the reader is 
referred to the studies mentioned above. The emission spectrum 
predicted by the model can be approximated as 

+ 5& [(! - y) ( mode1 ' U **J + 7 Mmodel, U mia , t/ max )] (3) 

where the first term accounts for the residual infrared emission 
coming from the stars and is given by the product of Q star (the 
solid angle subtended by the stars, and B v (T sVtt ), the blackbody 
emissivity with a fixed temperature r star = 5000 K, which was 
found by Smith et al. (2007) to provide a good description of 
the stellar continuum for A > 5 pirn. D is the distance of M3 1 
(we assumed 778 kpc), jy (model, U m i n ) is the dust power ra- 
diated per unit frequency per H nucleon when the dust is ex- 
posed to a single radiation field of intensity U m i n . The factor y 



(0 < y < 1) is introduced to parametrize the dust heating ef- 
fects of a power-law distribution of starlight intensities, as de- 
scribed in Draine & Li (2007). Actually, j' v (model, £/ m ; n , i/ max ) 
gives the dust power radiated per H nucleon by dust exposed to 
a power law distribution of starlight intensities (oc U~ 2 ) com- 
prised between U m \ n < U < U m!ix . The fraction y (0 < y < 1) of 
dust emission associated with this intensity field represents dust 
emission close to OB associations and/or photodissociation re- 
gions, where the intensity U > 1/^; in such a way this approach 
allows to handle dust temperature variations, in a smooth and 
convenient way. Following Draine et al. (2007) we fixed in our 
calculation i/ max = 10 6 . Finally Mh is the total mass of hydrogen 
(providing that the models give the emission per H nucleon). The 
mass of the dust, Md us t, can be obtained from Md us t = Mu (-^p), 
assuming the dust to gas ratio of the models. 
Each spectrum is thus completely characterized by 5 free pa- 
rameters: the solid angle subtended by the stars Q stal , the kind 
of emission model (amount of PAHs), the minimum intensity 
U m i n of the starlight radiation field, the fraction y of the total 
dust mass heated by the power-law distribution of starlight in- 
tensities, and the total dust (hydrogen) mass Mh- The best fit is 
obtained minimizing the quantity: 

2_ 1 y (Fv,obs,b — (^v.model)) 2 ... 
b obs.fr model 

In Table 2 we collected all the infrared integrated flux mea- 
surements of M3 1 performed so far, obtained with different in- 
struments: IRAC (Barmby et al. 2006), MIPS (Gordon et al. 
2006), COBE (Odenwald et al. 1998), IRAS (Rice et al. 1988), 
MSX (Kraemer et al. 2002) and ISO (Haas et al. 1998). In the 
fourth column we show our own measurements derived from 
IRAC and MIPS observations. In general all these measurements 
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are in good agreement, and define the spectral energy distribu- 
tion of M3 1 from ~ 1 /urn up to ~ 250 //m, as shown in Fig. 3. 
The best fit parameters we obtained were: PAHs = 4.6%, 
U min = 0.4, y = and M dust = 7.6x 10 7 M Q , Q star = 3.1 x 10~ 16 sr 
and the^ 2 of the best model was^ 2 = 0.55. 

In the following we discuss the uncertainties on the derived 
parameters and the implications of the results. 

3.1. The uncertainty in the dust mass estimate 

The mass estimate obtained in the above paragraph must be con- 
sidered rather uncertain. This is due to the fact that, as shown 
in Fig. 3, the peak wavelength of M31's IR spectrum is close to 
160yum. Thus, the vast majority of our analyzed observations 
is found at shorter waveleghts, therefore colder mass compo- 
nents can't be reliably constrained by the present observations. 
The only measurement that gives an appreciable constraint is 
the COBE 248 fim observation. Otherwise, the cold-mass esti- 
mate critically depends on the temperature of the dust inferred 
from the infrared spectrum. For a given mean radiation field in- 
tensity U at which the dust is exposed the steady-state temper- 
ature of the dust T °c t/ 1/6 and the emission peak A cc U~ 1/b 
(for grains with dimension > 30A, Draine & Li 2007), that 
is dT/T cc 1/6 dU/U. By energy balance considerations U is 
then directly proportional to the emissivity per H nucleon of the 
dust grains tabulated in the models and consequently the dust 
mass Mdust derived comparing the model with the observations 
is Mdust K U l , or fiMdust/Mdust <* -dU/U. This makes the dust 
temperature estimate relatively robust against uncertainties in U, 
i.e. a variation of 100% in U (and in Md us t) produces a vari- 
ation of only ~ 17% in dust temperature and in emission peak 
wavelength, to which our measurements are sensitive. Moreover, 
dU/U can be quite large even for small changes of U if U is also 
small, as it seems to be the case for M3 1 . In order to perform a 
more quantitative analysis, in Fig. 4 we show the^ 2 of the best 
model obtained fixing the value of t/ m ; n , against the correspon- 
dent total hydrogen mass M H predicted by the model. The black 
points (continuous line) in that figure show the result of the fit 
when all measurements in Tab. 2 are considered (with the ex- 
clusion of ISO and MSX measurements which were not used 
here), the red points (dotted line) when neglecting the COBE 
248 fim, and the blue (dashed line) points when neglecting the 
MIPS 160 fim and the COBE 148 fim and 248 fim measurements. 
From that is clear that the cold mass can't be accurately esti- 
mated. All the curves show a shallow minimum in x 2 but the 
range of mass values at which x 2 approaches the absolute mini- 
mum is large even considering the fit with all the present obser- 
vations. Moreover, once FIR measurements are excluded the^ 2 
curves are even flatter increasing the uncertainty in the mass val- 
ues. Nevertheless it has to be noted that all^- 2 curves are highly 
asymmetric and imply that large mass values are much less con- 
strained than low mass as explained below. 

In Fig. 5, we present the difference between the observa- 
tions and the values predicted by two different models: (i) the 
best fitting model (which has t/ m ; n = 0.4, black filled points); 
(ii) the best model obtained imposing t/ mm = 2 (red open cir- 
cles). All the measurements reported in Tab. 2 are considered 
in the fit, thus the x 2 of these two models are shown in Fig. 4 
for the correspondent values of t/ m ; n as given by the points 
connected by the continuous black line. In Fig. 5, only mea- 
surements with A > 50 fim are considered because this is the 
spectral region where the major differences can be noticed. A 
model should be considered consistent with the observations if 



all the observation minus model differences (A) are consistent 
with zero, within the observational errors. Focusing the atten- 
tion on the Far Infrared measurements (FIR) of Spitzer 160//m 
and COBE 148 fim, 248 fim in Fig. 5, we obtained for the best 
fitting model A 148/Jm = 0.55 cr 148/jm , A 160flm = 0.24 cr 160/jm and 
A248//111 = -0.06 a"248 J um, whereas for the model with t/ min = 2 
we obtained Ai 48jUm = 1.74 cr 148/jm , A 160/im = 1.57 cr 160/jm and 
A24 8//m = 2.45 cr 248/jm . Thus in this spectral range the best fit- 
ting model with t/ m ; n appears consistent with the observations 
at < 0.6 cr, whereas to reconcile observations and model pre- 
dictions once t/ m i n = 2 we have to admit differences < 2.5 cr, 
and in any case > 1 .5 cr. For smaller wavelengths the two mod- 
els appear almost equivalent in reproducing the observations, al- 
though a slightly worse result is obtained for t/min = 2. An im- 
portant feature of Fig. 5 is that for A > 140 /im the flux model 
values for t/ m ; n = 2 are lower than the values of the best fitting 
model, whereas for A < 140 fim they are systematically larger. 
This is what we expect as the peak emission of the model with 
t/min = 2 is shifted towards smaller wavelengths with respect to 
the best fitting model, as a consequence of the larger radiation 
field that is heating the dust. Because of model predictions in the 
FIR for the case of t/ m i n = 2 are already at the limit of predict- 
ing an acceptable wavelength distribution, larger radiation fields 
can be discarded. We conclude that on the basis of the diffuse 
dust emission models of Draine & Li (2007) and the analysis of 
the FIR spectrum of M3 1 the mean radiation field 3 that is heat- 
ing the dust in M31 is typically U < 2. As the radiation field 
is inversely proportional to the dust mass estimate we can also 
conclude that the mass of the dust in M31 is > 1.1 X 1O 7 M , 
which is the value obtained for the best model with t/ mm = 2, 
assuming a dust-to-gas ratio equal to 1%. 

The most recent and accurate estimate of the total mass 
of atomic hydrogen (HI) in M31 is that one of Braun et 
al. (2009), and is equal to M m = 7.33 x 10 9 M o . Nieten et 
al. (2006) provided a mass of molecular hydrogen Mh, equal 
to Mh 2 = 3.6 x 10 8 M o (within a radius of 18kpc). These 
recent estimates combine to a total neutral hydrogen mass of 
M ra + Mh 2 = (7.69 . . . 8.05) X 10 9 M Q once considering an un- 
certainty range for the CO conversion factor between X m = 
(2. . .4) x 10 20 molcm _2 (Kkmsec _1 ) _1 . The mass estimate we 
obtained from the best fitting model (considering all the ob- 
servations in Tab. 2, Mh = 7.6 x 10 9 M o ) is thus respectively 
1% and 5% smaller with respect to the neutral hydrogen mass 
estimates at the extremes of the uncertainty range mentioned 
above. Thus despite the large uncertainty we dicussed previ- 
ously, the best fitting model mass estimate appears consistent 
with HI and CO measurements, as shown in Fig. 4. Earlier esti- 
mates of the dust mass (M<j ust ) in M3 1 provided values equal to 
M dust = 3.8 xlO 7 M (Haas etal. 1998) and M dus t = 1.3xl0 7 M o 
(Schmidtobreick et al. 2000) which are fully consistent with our 
lower limit estimate but, for what have been said so far, they 
cannot be considered more accurate than our best model value 
estimate. 

Draine et al. (2007) analyzing a sample of 17 galaxies from 
the SINGS-SCUBA sample have shown that in all these galax- 
ies the value of the mean radiation field derived from the best 
fitting model was > 2, with a median value of 4.3. On that basis 
they proposed a restricted fitting procedure for the case in which 
submillimetric observations are not available which would im- 



3 In principle the value of f/ m j n is not coincident with the definition 
of the mean radiation field given by Draine & Li (2007), although these 
two values are essentially equivalent once y ~ as seems the case for 
M31. 
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Fig. 3. Infrared spectral energy distribution of M31. Color symbols are measurements obtained with the different instruments indicated in the leg- 
end. The continuous line denotes the best fitting model of Draine et al. (2007). Open circles represent the model predicted fluxes after convolution 
with the instrumental response function, and these are the points that are compared with observations. 



ply to set t/ m in > 0.7 during the fit in order to avoid an overes- 
timate of the dust mass. Nevertheless as stated by the same au- 
thors this procedure underestimates the dust mass if the radiation 
field heating the dust is weak. The sample of SINGS-SCUBA 
galaxies could be biased towards increased associated star for- 
mation as recognized by the same authors. As we have demon- 
strated above, for the case of M3 1 models with U m i n > 2 are 
not in agreement with the FIR measurements. Moreover models 
with a strong UV radiation field imply gas masses far below the 
range given by neutral hydrogen mass measurements, whereas 
our best fitting model (C/ m ; n = 0.4) comes very close to the ex- 
pected value. This additional constraint rules out models with 
high U m \ n and proves that M3 1 is a galaxy where the mean UV 
radiation field heating the dust is very weak. In Sect. 4 we will 
provide an independent demonstration of this fact and discuss its 
implications. 

In Fig. 6 (upper panels) we show the correlation between 
the logarithm of the total IR (TIR) emission and the logarithm 
of the 24/im (160yum) emission. The TIR was estimated from 



the S pitzer maps of M3 1 using the formula given by Draine & 
Li (2007): 

/tir = 0.95v/ y (8/zm)+1.15v/ v (24/im)-i-v7 y (70 i um)+v7 v (160 J um)(5) 

The lg(160yum) - lg(TIR) correlation is less scattered 
(RMS = 0.03) with respect to the lg(24jt/m) - lg(TIR) corre- 
lation (RMS = 0.06). We fitted these data with a simple linear 
model obtaining: 



Ig(TIR) = 0.89(0.02) lg(vF" y s [24fim]) - 0.1(0.2) 
lg(TIR) = 1.1(0.01) ^^''[160^ m])- 1.3(0.1) 



(6) 
(7) 



Although, as shown below, the correlation between the 24 yum 
(160 fim) and the TIR emission is not strictly linear, these results 
(in particular the calibration relation between the lg(24 fim) and 
the lg(TIR) variables) will be useful later in Sect. 5 and suffi- 
ciently accurate for our purposes. 

In Fig. 6 (bottom panels) we show the 24yum/TIR 
(160 yum/ TIR) ratios against the TIR intensity. As the TIR 
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Fig. 4. Minimum x 2 against Mu derived from the model fits, for each 
fixed value of (/ mill indicated by the labels. The black points (continuous 
line) represent the fitting results when considering all the measurements 
in Tab. 2 (with the exclusion of ISO and MSX measurements), the red 
points (dotted line) when neglecting the COBE 248 yum measurement, 
and the blue points (dashed line) when neglecting the MIPS 160yum 
and the COBE 148 yum, 248 yum measurements. The vertical dotted lines 
indicate the range of expected total hydrogen mass values estimated 
from HI and CO maps. Note that the points in this figure are connected 
by lines for clarity, but the^ 2 curve is actually more complex, showing 
local minima for each value of U mm . 



increases the 24yum/TIR ratio increases, but remains confined 
typically between 5%-7%. Conversely, the 160jum/TIR de- 
creases, changing from ~ 65% to ~ 50%. This demonstrates 
that overall the infrared spectrum of M3 1 is dominated by the 
cold dust emission in the FIR, although the contribute of hotter 
dust components tends to increase at larger TIR emissions. 

3.2. PAHs abundance, y 

Fig. 7, shows the ratio R$ of the 8 fim flux (normalized to the total 
infrared emission) against the ratio ^71 of 71 jjm flux to 160 /im 
flux. As proposed by Draine et al. (2007) this diagram can be 
used to estimate the abundance of PAHs. We used it to illustrate 
the uncertainty of our fit comparing the result obtained from the 
total integrated emission of M3 1 to that one obtained in the sin- 
gle regions of our maps, regridded at the 160/im resolution and 
using exclusively the regions with S /N > 3. Uncertainties were 
calculated as explained in Sect. 2. This selection restricted us 
to the study of the ring-spiral structure. Essentially PAHs are 
stochastically heated by single photons thus in order to derive 
their abundance one can compare the observations with theo- 
retical models with y — (Draine et al. 2007, where y is the 
parameter that regulates the power-law intensity radiation field 
in Eq. 3). In Fig. 7, we show the theoretical predictions as func- 
tions of the different amount of PAHs and of the minimum ra- 
diation field intensity U^ n . The results imply that models with 
large abundances (> 3.19%) and rather low t/ m j n are favoured 
when comparing to the data. In fact the measurement associated 
with the integrated infrared emission (blue point) is not at the 
center of the distribution of points coming from the analysis of 
the pixel-to-pixel analysis, and in particular is biased towards 
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ting model (U mm = 0.4), whereas red open circles considering the best 
model among those ones with U mm = 2 fixed. Errorbars are the uncer- 
tainties associated with each measurement as given in Tab 2. 

smaller values of Rg. As our analysis is restricted to the spiral- 
ring structure, this result suggests that in the inter-arm regions 
and/or toward the bulge the abundance of PAHs could be smaller 
than what we derived, but an accurate analysis of the distribution 
of these particles across the whole disk of M31 is beyond the 
purpose of this work. 

In Fig. 8 we show the ratio R24 of the 24 yum flux (normalized 
to the total infrared emission) against the ratio Rji,a similar dia- 
gram of Fig. 7, but for the 24 yum map, which is better diagnostic 
to the radiation field intensity (Draine et al. 2007). We consid- 
ered in this case only models with PAHs abundances > 3.19% 
since for the analyzed regions these are the models that better 
agree with observations (Fig. 7). For the bulk of the points, and 
for the total integrated flux measurement y can be in general 
considered < 0.05. Using the equation: 

f PDR = 1.05(*24 - 0.14*8 - 0.035) 75 (8) 

taken from Draine & Li (2007), and the observed ratios R& and 
R24 defined above, we determined the fraction of the total dust 
luminosity that is radiated by dust grains in regions with U > 10 2 
(/pdr) to be lower than ~ 4%. 

4. Which sources are responsible for the dust 
heating in M31? 

The analysis of the infrared spectrum presented in Sec. 3 allows 
to conclude that the infrared spectrum of M3 1 can be explained 
by diffuse dust emission. Nevertheless, on the basis of that, we 
cannot understand for which stellar populations are responsible 
for the dust heating process. In particular we want to know where 
the dust in M3 1 is predominantly heated by young populations 
( Age < 1 Gyr) or where older populations ( Age > 1 Gyr) play 
an important role in the dust-heating process as well. As shown 
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below the answer to this question is of crucial impact for the 
deduction of several important parameters, e.g. attenuation and 
SFR, and thus has to be addressed as a central issue. 

In the following we used the observational maps of M3 1 re- 
sampled at 40"/px as described in Sec. 2, and considered all the 
regions in which the S /N > 1 in all the maps. 

Also it should be noted that the results we obtained are based 
on the theoretical models of Kong et al. (2004) and Cortese et 
al. (2008). It is worth saying that alternative interpretations could 
exist of the observational data we present. In particular models 
considering modifications of the canonical attenuation laws have 
been proposed (e.g. Inoue 2005, 2006 see later in this Section). 

The model of Kong et al. (2004) is based on the Bruzual & 
Chariot (2003) population synthesis code, it assumes an expo- 
nential star formation history (SFH) and a power-law absorp- 
tion law, distinguishing between young (< 10 7 yr) stars embed- 



ded in their birth clouds and older stars migrated into the ISM. 
The model of Cortese et al. (2008) is also based on the Bruzual 
& Chariot (2003) population syntesis code, but is assumes an 
LMC attenuation law, and that both stars and dust were homo- 
geneously distributed in a plane parallel (sandwich) geometry. 
They also adopted a Star Formation History 'a la Sandage' in 
the formalism of Gavazzi et al. (2002): 

SFR{UT) = ±exp[-^ (9) 

where SFR is the Star Formation Rate per unit mass, t the age (in 
Gyr) of the galaxy (assuming t=13 Gyr at present epoch), and r 
is the time (in Gyr) at which the star formation rate reaches the 
highest value over the whole galaxy history (note that in this no- 
tation the present age of the stellar populations born at time t = r 
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is given by age =13 Gyr - t). As stated by the same authors, t 
should be considered as a proxy for the shape of the SED rather 
than being used to derive an exact indication of the age of the 
underlying stellar populations. We calculated for each value of 
r, the time t* where the SFR reaches half of its maximum value 
(for t > t) over the whole galaxy history. Also the age of the 
stellar populations correspondent to t = t* must not be consid- 
ered in a strict manner but it can be considered as an indication 
of the age of the youngest stellar populations that are signifi- 
cantly contributing to the SED, given that the contribute of next 
generations of stars becomes increasingly less important in or- 
der to explain the observed SED. For example, for t = 8 Gyr we 
obtain t* ~ 13.5 Gyr. Assuming for the present epoch t=13 Gyr 
as in Cortese et al. (2008), we can expect that young populations 
( age < 1 Gyr) contributed significantly to the observed SED. 
For t = 5 Gyr we obtain instead t* ~ 9.6 Gyr which implies 
that stellar populations contributing to the SED should have an 
age > 3 Gyr. 

It should be also remind that the predictions of these models 
should not be considered reliable in regions where old stars are 
likely to dominate the UV flux (mostly in the bulge region for 
M31). 

At first we investigated the TIR/FUV vs. (3 (the UV slope) 
relationship for the analyzed regions as shown in Fig. 9. The 
solid line represents the fit obtained by Kong et al. (2004) for 
the 50 starburst galaxies in their sample, which represents the 
so called IRX - /3 relationship for starburst galaxies (Meurer et 
al. 1999). We adopted the same definition of Kong et al. (2004) 
to calculate /?: 



(10) 



board of GALEX and /fuV'/nuv ^ the mean flux densities 
(per unit wavelength) through these filters. 

The small black points are the values of the TIR/FUV and 
/3 obtained in the analyzed regions of M3 1 . A positive corre- 
lation is visible, but it does not follow the relation found in 
starburst galaxies. As can be expected naively, going towards 
negative slopes and bluer colors our observations imply simi- 
lar TIR/FUV and [} values with respect to the starburst galaxies 
relationship whereas redder regions behave differently showing 
typically lower TIR/FUV for equal /3. This is in line with Kong 
et al. (2004), who showed the results coming from the analysis 
of their sample of normal star-forming galaxies as well and ob- 
tained that these objects typically have lower values of TIR/FUV 
for fixed UV-color with respect to what predicted by the (IRX-/3) 
relation, with a large scatter of values. Other authors have re- 
cently questioned the application of the (IRX - /3) relationship 
for normal star-forming galaxies, suggesting that such systems 
suffer from lower dust attenuations with respect to what could 
be inferred from the (IRX - P) relation (Salim et al. 2007). 

Following the notation of Kong et al. (2004) we plot in Fig. 9 
the lines of constant ratio of present to past-averaged star forma- 
tion rate assuming an an exponentially declining star formation 
history (b parameter, blue short-dashed lines) and of constant 
attenuation (Auv, red long-dashed lines). 

This comparison suggests that in the M31 regions popula- 
tion gradients are present but there are similar amounts of dust 
attenuation, being our observational points elongated along the 
direction indicated by the bands of equal attenuation. 

As color variations seem due to a spread of ages, it is more 
convenient to adopt a larger color baseline like the (FUV-/sdss) 4 
color. Cortese et al. (2008) provided a detailed discussion of the 
age dependence of the TIR/FUV ratio, colors and attenuations. 



where Afuv - 1516 A and A^uv - 2267 A are the effective 
wavelengths of the far-ultraviolet and near-ultraviolet filters on 



4 We have decided to use here the (FUV - z'sdss) color as the 2MAS S 
images of M3 1 obtained during the 2MAS S Extended S urvey did not 
appear sufficiently deep for our purposes. 
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Fig. 8. Ratio of the 24 jim flux normalized to FIR emission against 7 1 fim to 1 60 yum flux ratio Rj 1 . The blue point is the ratio correspondent to 
the total integrated flux measurement of M31. Each horizontal red line shows the ratio values predicted by the models of Draine & Li (2007) for 
different values of the minimum radiation field intensity U mm and of the y parameter (see text). The three panels show the results for models with 
different abundances of PAHs as indicated by the labels. 



In accordance with their results, for (FUV - /sdss) ~ 4.3 the 
energy absorbed by the dust at A < 4000 A is approximately 
equal to the energy absorbed at A > 4000 A (as inferred from 
their Fig. 2 and Table 1). This corresponds to t* ^8.5 Gyr as 
defined above, thus the stellar populations that are contributing 
to the SED have an age> 4.5 Gyr. 

In Fig. 10 we show the TIR/FUV ratio against the (FUV - 
'sdss) color, for regions with S/N > 1 across the lOkpc ring 
(8kpc < r < 12kpc). The positive correlation between these 
variables is still clearly visible also in this diagram. Moreover 
the color range 3.5 < (FUV - 'sdss) < 7 demonstrates that in a 



large number of regions old populations dominate the dust heat- 
ing process. As shown in Fig. 10 (lower left panel), the mean 
Apuv in the analyzed regions of M3 1 derived with this method 
is Apuv = 2.1 ± 0.4 mag where the error is the standard devia- 
tion. If the starburst scenario is considered, we obtain the result 
shown in the lower right panel of Fig. 10, where the mean atten- 
uation is Apuv = 4.2 ±1.0 mag. Neglecting the dust heating age 
dependence would imply a mean attenuation 100% larger and 
a scatter 150% larger, since in this case the large color changes 
would totally be attributed to the reddening. 



M. Montalto, S. Seitz, A. Riffeser, U. Hopp, C.-H. Lee, R. Schonrich: Dust properties in M31.I. 



11 



3 



> 2 

D 

K 

i 1 

&J0 







A uv =4.5 








2 







Fig. 9. (TIR/FUV) logarithmic ratio against /3, the UV-slope, for each region with S/N > 1 analyzed in the 40"/px resolution maps. The thick 
continuous line denotes the best fit of Kong et al. (2004) for their sample of starburst galaxies. The long dashed red lines and the short dashed blue 
lines are constant attenuation and constant ratio of present to past-averaged star formation rate lines as given by Kong et al. (2004). 



In Fig. 10, we overplotted two lines of constant attenuation 
correspondent to Apuv = 2 and Apuv = 3 as indicated by the 
labels. The shape of these lines reflects the age dependence of 
the dust heating as modeled by Cortese et al. (2008). Actually 
for blue (FUV - /sdss) < 4) or red colors (FUV - /sdss) > 7 
the attenuation lines tend to be parallel to the color axis, and 
thus independent from the color. This is due to the fact that 
in these color ranges the intrinsic spectral energy distribution 
of the populations that are heating the dust is not expected to 
change appreciably with age to the very young (blue) or very 
old (red) side. On the other hand in the range comprised be- 
tween 4 < (FUV - /sdss) < 7, which is actually where most of 
our observations fall, age variations contribute significantly to 
color changes, and older populations appear intrinsically redder 
thus shifting the TIR/FUV ratio towards larger values for equal 
amounts of attenuation with respect to bluer regions. In the fol- 
lowing we tried to disentangle the effects of the age and of the 
reddening on the observed colors and TIR/FUV ratios. 

At first we converted the colors to the values of t (the time 
expressed in Gyr at which the star formation rate reaches the 



highest value over the whole galaxy history, assuming 13 Gyr at 
present epoch) by means of: 



lg(r) = -0.073 (FUV - /sdss) + 0.96 



(ID 



taken from Table 2 of Cortese et al. (2008). The colors in 
Fig. 10 are codified in function of the value of t. We obtained 
that 83% of the analyzed regions across the spiral-ring have 
(FUV - /sdss) > 4.3 and thus values of r < 4.4 Gyr. Hence 
according to the results of Cortese et al. (2008) in 83% of the 
analyzed regions the dust absorbs more than 50% of the energy 
at A > 4000 A. As stated above this implies that in these re- 
gions stellar populations responsible for the dust heating should 
be at least a few Gyr old. This could provide a good agreement 
to the low mean intensity of the radiation field we found from 
the analysis of the infrared spectrum of M31, although a small 
value of U could be obtained also by young stars with a low 
SFR. It is important to remind that these results do not imply 
that younger populations are not present at all, but rather that 
their contribute is not dominant in these regions. If, on the con- 
trary one tries to explain the observed colors assuming that only 
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Fig. 10. Upper panel: (TIR/FUV) logarithmic ratio against (FUV - i'sdss) color for regions with S/N > 1 across the 10 kpc ring (8kpc < r < 
Ylkpc). The colors denote the different values of t (the time expressed in Gyr at which the star formation rate reaches the highest value over the 
whole galaxy history, assuming 13 Gyr at present epoch) obtained for each region following the recipes of Cortese et al. (2008). The two solid 
lines are two lines of constant attenuation. The errorbar on the top-left corner shows the mean observational errors. Lower panels: histograms of 
ultraviolet attenuations derived following the recipes of Cortese et al. (2008, left panel) and assuming the starburst scenario (right panel). Mean 
values and 1 a uncertainties are denoted by solid and dashed vertical lines and are A FUV = 2.1+ 0.4 mag and A FUV = 4.2 ± 1.0 mag for left and 
right panel respectively. 
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young populations (age< 1 Gyr) are heating the dust, one should 
also admit much larger attenuations for normal galaxies, as de- 
mostrated above and for the same distribution of star and dust 
as in Cortese et al. (2008). Otherwise it is clear that in order to 
confirm these results more accurate models are needed. Some 
other considerations regarding the limitation of the Cortese et 
al. (2008) approach are presented shortly below and in the next 
Section. 

From another point of view we show in Fig. 1 1 the same 
diagram of Fig. 10 separating regions of low and high TIR 
emission as indicated by the labels and the colors. We consid- 
ered 4 equal bins of TIR emission comprised between 10~ 11 - 
10~ 10 ergcirT 2 sec -1 . From Fig. 11 it appears that moving to- 
wards redder colors along the regions of constant TIR emission, 
the TIR/FUV ratio increases. As the TIR emission is constant, 
the FUV emission must decrease. Nevertheless, if the lower 
FUV emission was due to larger attenuations, the TIR emis- 
sion should have increased accordingly, which is not observed. 
Thus, the color change must be due to an intrinsic lower emis- 
sion in the FUV band, that could be indicative of the presence 
of a population gradient. Strictly speaking this interpretation is 
valid for an homogeneous layer of stars embedded in an homo- 
geneous optically thin layer of dust. In a real galaxy dust and 
stars (especially recently formed stars) are likely to have more 
clumpy and inhomogeneous distributions. Otherwise, studying 
two of the largest dust clouds listed in the M31 Atlas (Hodge 
1980), Hodge & Kennicutt (1982) derived average extinctions 
of Ab = 0.35 mag (with a maximum value of Ab = 0.60 mag) 
and A B = 0.43 mag (maximum A B = 0.76 mag) for D307 and 
D441 respectively, indicating that the majority of all the other 
clouds should have smaller amounts of extinctions. Other studies 
confirm the idea that the interstellar medium in M3 1 is generally 
optically thin. Fan et al. (2007) derived reddening values towards 
443 M3 1 's catalogued globular clusters obtaining that more than 
half of them are affected by a reddening E(B - V) <~ 0.2 mag 
with an average value of E(B - V) — 0.28^ ^ (see also Barmby 
et al. 2000). Regarding the homogeneity of the dust distribution, 
Inoue (2005, 2006) investigating the attenuation law resulting 
from clumpy spatial distributions of dust and (young) stars, pro- 
posed different explanations for the redder UV colors of normal 
galaxies with respect to starburst galaxies for fixed FIR to FUV 
ratios, which are not based on population gradients. These al- 
ternative models can be summarized as: (i) models where the 
attenuation law has a steeper dependence on wavelength than 
the canonical attenuation laws as a result of an age-selective at- 
tenuation of young stars in clumpy structures with respect to 
old stars; (ii) models where the attenuation law has no atten- 
uation bump at 2 175 A; (iii) models with a bump at 2 175 A, but 
with smaller albedo for shorter wavelengths (except for the bump 
range). While we did not try to apply these models to our data, 
we observe that at least the first two scenarious do not seem to fit 
the case of M31. Bianchi et al. (1996) derived that the M31 ex- 
tinction curve is very similar to the average Galactic extinction 
law and a possible reduction of the 2 175 A bump is significant 
only at the 1 <j level. Similar results were reached also in previous 
studies (e.g. Walterbos & Kennicutt 1988). As for the third class 
of models further studies of the wavelength dependence of the 
albedo are necessary. We also observe that these considerations 
do not imply that clumpy structures with young stars embed- 
ded do not exist in M31, but it could be that most of the dust in 
M31 is not located in these structures. From counts of dark neb- 
ulae across the disk of M3 1 Hodge & Kennicutt (1982) obtained 
that the major dust lanes visible in the optical should account for 
~ 15% of the total dust mass content. Recent results by Nieten et 



al. (2007) show a good correlation between the most prominent 
dust lanes and dense molecular clouds traced by strong CO emis- 
sion lines but concluded that molecular gas in M31 is only ~ 7% 
of the total neutral gas content in M31, and that dust appears 
correlated also with atomic gas in M31 which is distributed in 
more extended regions out of the densest clouds in the spiral-ring 
structure. Moreover in Sec. 3.1 we derived that only a fraction 
< 7% of the TIR comes from the hot dust component emission 
at 24/im associated with young star forming regions. Finally, it 
is unlikely that a large fraction of the M3 1 's dust mass is in very 
cold dark clouds (T < 16 K), as the total M31 dust mass we 
derived in Sec. 3.1 from dust diffuse emission models was close 
(l%-5%) to the estimate based on neutral gas measurements. We 
conclude that the interpretation of Fig. 1 1 in terms of population 
gradients must be considered carefully and in need of further in- 
vestigation, but it appears reasonable for regions dominated by 
diffuse dust emission away from large star-forming complexes 
(in this regard remember the criteria we applied in Sec. 2 against 
the selection of HII regions). In fact the attenuation along bands 
of equal TIR emission appears to decrease towards redder colors 
accordingly to the predictions of Cortese et al. (2008). This could 
be also in agreement with the idea of a population gradient, as 
older populations would tend to be intrinsically less attenuated 
than younger ones, because of their intrinsically redder spectral 
energy distribution, at least in the region where the age effect is 
maximal 4 < (FUV - /sdss) < 7. We note nevertheless that the 
observed color could be also biased towards redder values es- 
pecially in the bulge region where the geometry of the dust and 
star is clearly deviating from the simplyfied sandwich model of 
Cortese et al. (2008). We will discuss such potential bias effect 
in more details in the next Section. Another intresting consid- 
eration regarding the diagram of Fig. 11, is that moving verti- 
cally towards larger TIR/FUV ratios the TIR emission increases 
(whereas the FUV emission remains almost constant), which 
should imply larger attenuation values because of the large ab- 
sorbed energy for a fixed observed FUV energy, which is what 
the model predicts. We have thus the following interpretation of 
the diagram in Fig. 1 1 . While the large color variation of the an- 
alyzed regions seems to be due to the gradient of populations 
that are heating the dust, the spread of the TIR/FUV values for 
a fixed color seems instead due to the differential reddening of 
regions with similar underlying stellar populations. 

We further investigated the spatial distribution of the A FUV 
across the disk of M3 1 . We considered nine 2 kpc-wide radial 
bins from the center of M31 up to 18kpc as shown in Fig. 2, 
and imposed a S /N > 1. In Fig. 12 we present the logarithmic 
TIR/FUV ratios against the (FUV - /sdss) colors obtained in 
each radial bin. In general it appears that these observable quan- 
tities are positively correlated across all the disk of M3 1 . Only in 
the innermost annulus it was not possible to verify the positive 
correlation of these variables because of the few regions ana- 
lyzed. We fit the lg(TIR/FUV) ratios against to (FUV - / SD ss) 
colors in each radial bin with a linear model: lg(TIR/FUV) = 
a (FUV - /sdss) + b . We considered separately the uncertainties 
along the x and y axis and took the mean value of the derived pa- 
rameters and their semi-difference as estimates of the best model 
parameters (a and b) and uncertainties. In Tab. 3 we reported the 
result of the fit along with other useful quantities. 

As shown also in Fig. 13 the mean attenuation in each radial 
bin varies reaching a maximum near the lOkpc ring, while it 
seems to decrease faster towards the inner regions of M3 1 than 
in the outer regions. 

These results are also shown as the color maps of Fig. 14. 
While as dicussed above the model adopted here could be less 
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Fig. 11. TIR/FUV logarithmic ratios vs. (FUV - /sdss) colors for regions with S/N > 1 across the lOkpc ring (Skpc < r < \2kpc). Different colors 
and labels denote regions of constant TIR emission. Red, blue, black and magenta lines indicate the best fitting linear regression models of the 
data in each TIR emission interval. Green curves are two constant attenuation lines accordingly to the recipes of Cortese et al. (2008) as shown in 
Fig. 10. 



reliable in the inner part of the galaxy due to deviations of the 
dust and star distribution from the assumed geometry, the trend 
and the values of the mean attenuations derived with this method 
are in agreement with previous studies as discussed in the next 
Section. 

The existence of a correlation between the TIR/FUV ratios 
and the colors that we have found needs a final remark. Actually 
it is important to keep in mind that the theoretical calculations 
of Kong et al. (2004) and Cortese et al. (2008) do not neces- 
sarily imply the existence of any correlation between these ob- 
servables. Nonetheless we should remind that in this case we are 
studying regions inside one galaxy, and thus we can expect a 
more uniform and homogeneous behaviour. 

5. Discussion 

The models of Cortese et al. (2008) assume a simple plane paral- 
lel (sandwich) geometry for the dust and star distribution. Even 
if the fractional scale height between the dust and the stars is al- 



lowed to vary with wavelength no radial dependence is assumed. 
The method solely grounds on the relationship between two ob- 
served emission ratios and can by this be expected to be highly 
prone to light contamination from distant sources and a different 
geometry. Since, if the geometry is not plane parallel the radia- 
tion field experienced by the dust can be significantly different 
from what is naively measured by the outside observer, the ob- 
served color from which the age of the underlying population 
is estimated could be biased towards redder values in regions 
where older stellar populations are less embedded in the dust. 
This bias should be strongest in the inner regions of the galaxy 
as the result of the bulge light contamination. Regions with 
r < 8kpc in Fig. 12 reach redder colors and higher TIR/FUV 
ratios with respect to the outer regions. For r > 12kpc the color 
range remains almost constant although the TIR/FUV ratio de- 
creases. In presence of the above mentioned observational bias 
the value of the attenuation derived from the model could have 
been underestimated in the inner regions. Moreover in the inner- 
most radial bin (r < 2 kpc) old stars can contribute a significant 
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Fig. 12. Logarithmic TIR/FUV ratios against FUV - i'sdss colors for regions inside each annulus. Thin black solid and dashed lines are best fitting 
linear models of the data as explained in the text. Green curves are constant attenuation lines accordingly to the recipes of Cortese et al. (2008) as 
shown in Fig. 10. 



fraction of the GALEX-UV emission, and thus the model pre- 
diction should not be considered reliable in that region. 

To estimate how the integrated light along the line of sight 
could have affected our analysis and consequently biased the de- 
rived attenuations we compared the (FUV - !sdss) color along 
the major and minor axis as shown in Fig. 15. Each point in 
these figures is the average value of the colors (attenuations) in- 
side a 200 -wide (0.75 kpc) stripe centered on the major and 
minor-axis in correspondence of each one of the 2kpc annular 
regions shown in Fig. 2. The bias effect should be more evident 
along the minor than the major axis because of the inclination 
of M3 1 . At the same galactocentric radius of the disc the line 
of sight passes through more central regions of the bulge along 
the minor axis than along the major axis. In case the described 
bias affects the derived parameters it has thus to be expected that 
at the same galactocentric radius the adopted method will pre- 
dict lower attenuations along the minor axis than on the major 
axis. We calculated the mean color and attenuation difference 
inside the lOkpc ring excluding the region around +5 kpc on the 
near-side of the galaxy along the minor axis because of the pres- 



ence of the inner spiral arm with enhanced TIR emission and 
attenuation. In Fig. 15 (left panel) we see that indeed the color 
profile along the minor axis is systematically redder than along 
the major axis on average of ~ 0.5 mag. This implies that the 
attenuations appear to be systematically lower on the minor axis 
of A FU v ~ 0.2 mag (Fig. 15 right panel). 

Because also the major axis should be affected by such ob- 
servational bias the effect should be in principle larger than this. 
Considering the radial bin at 1 lkpc and that one at 3kpc and con- 
sidering the color shift toward redder values in the inner bin as 
totally due to the bias effect we estimated the average attenua- 
tion at 3kpc to be larger of ~ 0.84 mag (Fig. 16) than what de- 
rived from the observed colors. This value should be considered 
as an upper limit as the color shift should be in part attributed 
also to intrinsically redder and older populations than those in 
the ring where most of the star formation in M3 1 is taking place. 
Nevertheless even applying such a generous bias correction the 
ring would still stick out as a maximum of attenuation as shown 
in Fig. 16. On the other hand such an exercise points out the need 
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Fig. 13. Average FUV (A^ = 1516A) attenuations in 2kpc radial bins from the center of M31 to ~21kpc obtained by us (filled circles), Xu & Helou 
(1996, open circles). Open triangles represent the result of Hodge & Lee (1998) obtained analyzing 5 fields at different galctocentric distances. 



of a more accurate model to properly account for what happens 
in the innermost regions. 

Xu & Helou (1996) developed a dust heating/cooling model 
for M31 based on a radiative transfer code which assumes 
a sandwich geometry for stars and dust like in Cortese et 
al. (2008). They obtained that the mean optical depth t v viewed 
from the inclination angle of 77° increased with radius from 
t„ ~ 0.7 at r = 2 kpc, reaching a peak of t v = 1 .6 near 10 kpc and 
remaining quite flat out to 14 kpc. We used Eq. 2, 3, 4 of Cortese 
et al. (2008) to convert the t v values of Xu & Helou (1996) into 
Apuv magnitudes. The resulting values are Apuv = 1.23 mag 
at 2 kpc, and Apuv = 2.50 mag near 10 kpc in good agreement 
with the results shown in Tab. 3 and Fig. 13. The radial trend 
we derived actually peaks close to 10 kpc and appears to re- 
main quite flat in the outer regions and to decrease faster in 
the inner regions. While, considering the uncertainties in the 
Apuv estimate, these results are consistent with those ones of 
Xu & Helou (1996), some differences are visible. In the inner 
side of the ring our average attenuations tend to be lower (from 
~ (0.1 .. . 0.5) mag) with respect to the results of those authors. 
Despite the fact that the same sandwich geometry is adopted, the 



two approaches actually differ in various aspects. At first Xu & 
Helou (1996) did not admit the presence of population gradients. 
While this assumption was motivated by the consideration that 
the (V - R) color of M31 appeared rather constant in the disk 
of M31 (Walterbos & Kennicutt 1988) the results shown in the 
previous section, based on a much larger color baseline, indi- 
cate that indeed a population gradient is present across the disk 
of M3 1 . In absence of population gradients, as discussed in the 
previos section, one tends to overestimate the attenuation in the 
region where the age effect is maximal (4 < (FUV - z'sdss) < 7), 
and this independently on the observational biases discussed pre- 
viously (which should affect both methods in a similar way). It 
is interesting to note that in the inner region at 2kpc, and thus 
where most of the regions have (FUV - /sdss) > 7 the two mod- 
els appear in perfect agreement. On the other hand the model of 
Xu & Helou (1996) actually solves the radiative transfer prob- 
lem and fully accounts for the effects of the dust scattering which 
are neglected by Cortese et al. (2008). If the scatter is neglected 
for a fixed observed color we would have a smaller TIR/FUV 
ratio and thus we would predict a smaller attenuation. While 
both these effects could have contributed to the differences we 
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Fig. 14. From top to bottom: (FUV - z'sdss) map, logarithmic TIR/FUV map, and Apuv attenuation map. Color bars indicate the range of values 
spanned by the quantities using a linear scale. The first two maps are derived combining the FUV, ; S dss and all the infrared maps used in this 
work. The regions shown here have S /N > 1, FUV/TVf/V > 0.3 and 1 < IF^. < 10, and avoid bright sources detected in the ;'sdss band 
(black dots). The last map shows the attenuations A FUV derived using the recipes of Cortese et al. (2008) and the two above maps. The attenuation 
is maximum in the lOkpc ring where it reaches a mean value of 2.5 ± 0.7 mag, and declines more strongly towards the inner regions of the galaxy 
than towards the outer regions. 



see in Fig. 13 the general good agreement of these results telling 
us that these effects are not dramatically changing our conclu- 
sions. As we have shown in the previous sections old popula- 
tions are overall dominating the radiation field, so good results 
can be obtained if the color is fixed to the one of an old popula- 
tion as done in Xu & Helou (1996). The scattering should affect 
mostly young stars embedded in the dust, and thus mainly in the 
spiral-ring structure, but as the radiation field of these stars is 
never dominant the effect should be limited. Finally we should 
remind that we are comparing results obtained at different spa- 
tial resolutions (40" /px in our case and 2'/px in the case of Xu 
& Helou 1996) and that the accuracy of Spitzer, GALEX and 



SDSS maps is certainly larger than the measurements used by 
Xu & Helou (1996). In particular using IRAS images the cold 
dust emission and thus the TIR emission should be less accu- 
rately constrained than using Spitzer data as the longest wave- 
length measurement of IRAS is at lOOyum whereas for Spitzer 
is at 160 jt/m close to the FIR spectrum peak of M31 (Fig. 3). In 
Fig. 13 we show also the comparison of our results with those of 
Hodge & Lee (1998) who used two color diagrams to determine 
the average reddenings of 5 fields at different M3 1 galactocentric 
distances. We converted the reddening values of Hodge & Lee 
(1998) to the optical depths along the line of sight (t v ) by mul- 
tiplying a factor of 2x0.921 x 2.8 as in Xu & Helou (1996), then 
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Fig. 15. Left panel: (FUV - j'sdss) color profiles along the minor (red line) and major axis (black line). Right panel: average attenuations along 
the minor and major axis. The profiles are derived considering regions with S/N > 1. Each point is the average value of the colors (attenuations) 
inside a 200 -wide (0.75 kpc) stripe centered on the major and minor-axis in correspondence of each one of the 2kpc annular regions shown in 
Fig. 2 




Fig. 16. Upper panel: observed colors and TIR/FUV ratios for regions in the 3kpc radial bin (big red points) and in the 1 lkpc radial bin (small blue 
points). Lower panel: same as upper panel, but the observed colors in the inner radial bin are shifted by the average color difference with respect 
to the best fit linear model in the outer ring (blue line). 



we proceeded as before to convert the r v to Apuv- Although the 
radial trend always peaks at lOkpc a larger scatter with respect to 
the results of Xu & Helou (1996) is visible. This may be anyway 
accounted by the fact that both us and Xu & Helou (1996) are 
analyzing radial averages whereas the results of Hodge & Lee 
(1998) are based on 5 selected regions. 

The lOkpc ring is the locus where also gas density and 
star formation activity are maximum across the disk of M31 
(e.g. Braun et al. 2009; Nieten et al. 2006; Reddish 1962; 
Hodge 1979; Hodge & Lee 1988) thus is not surprising to ad- 



mit that also the dust density peaks at the ring as obtained also 
in previous studies. 

Another limitation of the approach we followed consists in 
the assumption that each one of the cells we analyzed is indepen- 
dent from the others, and thus the radiation field of each region 
is the only resposible for the dust heating in that region. As the 
mean free path of photons is typically larger than the dimension 
of the regions here anayzed ( 40 /px, ~ 670 pc/px), it is possible 
that photons coming from nearby bright regions could contribute 
significantly to dust heating. Nevertheless this effect is likely to 
have a limited impact on our conclusions. At first because dust 
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heated by stars in star forming regions is not the dominant source 
of emission in the FIR in M31, as given by the fact that the hot 
dust component emission at 24//m is overall < 7% of the TIR 
emission as shown in Sec. 3.1. Second of all the smooth distri- 
bution of the M3 1 surface brightness in the optical (A > 4000A 
from where most of the dust-heating radiation comes) suggests 
the hypothesis of homogeneous conditions among adjacent re- 
gions for which a local equilibrium between photon transmitted 
and receiveid in closeby regions can be assumed. 

Finally, using the correlation between the 24 //m and the TIR 
emission shown in Fig. 6 we derived the reddening map of M31 
at 6"/px (Fig. 17). We then applied the reddening correction to 
the observed SDSS map (Fig. 18). The most prominent dust 
lanes and absorption features are clearly addressed by the red- 
dening map, and the ring emission at lOkpc gains importance. 
Note that in Fig. 17, 18 we show the whole maps for clarity, but 
regions where optical, UV and infrared point-like sources are 
visible are not reliable. We discussed the issue of contaminants 
and how we corrected for them in Sec. 2. 

6. Conclusions 

In this paper we have explored some properties of the dust in 
M3 1 . The major results we reached are the following: (i) from 
the study of the infrared spectrum of M3 1 we have obtained that 
the mean intensity of the radiation field that is heating the dust 
is globally low (typically U < 2); (ii) the dust mass (Md us t) es- 
timate remains uncertain due to the lack of submillimetric ob- 
servations, but we have obtained that Md us t S 1.1 x 10 7 M©, 
the value given by the best fitting model being 7.6 x 10 7 M o 
in good agreement with what inferred from CO and HI obser- 
vations; (iii) the abundance of Polyciclic Aromatic Hydrocarbon 
(PAH) particles in M31 is high (> 3% M^t) across the spiral- 
ring structure of M3 1 ; (iv) we demonstrated the existence of a 
correlation between the observed TIR/FUV emission ratios and 
the color (FUV - /sdss) overall the spiral-ring structure of M3 1 ; 
(vi) this correlation is not in agreement with the IRX-/3 relation- 
ship of starburst galaxies, thus color changes are in general not 
driven by dust attenuations; (v) we found that according to the 
prescription of models which consider the age-dependent dust 
heating the observed correlation could be explained as the ev- 
idence of the presence of a population gradient and of a quite 
homogeneous attenuation of the analyzed regions; (vi) in 83% 
of the regions comprised between 8kpc < r < 12kpc the dust 
absorbs more than 50% of the energy at A > 4000 A, and it ap- 
pears mainly heated by populations a few Gyr old and this could 
provide a good interpretation to the low mean intensity of the 
radiation field we found from the independent analysis of the 
infrared spectrum; (vii) we determined that the mean attenua- 
tion on 2kpc-wide radial bins reaches the maximum value near 
lOkpc and decreases faster in the inner than in the outer regions 
of the galaxy, and that regions with larger TIR emission have 
also enhanced attenuations; (viii) finally we derived an attenu- 
ation map of M31 at 6 /px resolution (~ lOOpc/px along the 
plane of M31). 

Future contributes will investigate in more detail the depen- 
dence of the dust attenuation values from other parameters (e.g. 
color, metallicity), the star formation rate and the modellization 
of the stellar populations in M3 1 . 

The maps presented in Fig. 14 and the E(B-V) map of Fig. 17 
are avaiable via CDS. 
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Fig. 18. Upper panel: original SDSS image in the i-band. Lower panel: SDSS i-band image deredenned with our extinction map. For clarity all 
regions have been shown here. Regions with bright point-like sources and/or with occasional negative fluxes in some (UV/Optical/IR) maps (black 
regions) are not reliable. The colorbar shows intensity levels in Jansky in logarithmic units. Images have 6 "/px resolution. 
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Table 2. Measurements of M3 1 infrared integrated fluxes obtained using IRAC (Barmby et al. 2006), MIPS (Gordon et al. 2006), COBE (Odenwald 
et al. 1998), IRAS (Rice et al. 1988), MSX (Kraemer et al. 2002) and ISO (Haas et al. 1998) instruments. The last column shows our own 
measurements derived from IRAC and MIPS observations. We measured the integrated flux in an ellipse of semi-major axis 84' and semi-minor 
axis 18' centered on M31 with a position angle of PA = 38°. 



X" 


x a . 

mm 


max 


This work 


IRAC 


MIPS 


COBE 


IRAS 


MSX ISO 


/im 




fim 


Jy 


Jy 


Jy 


Jy 


Jy 


Jy Jy 


1.27 


1.1 


1.5 








534 ± 107 






2.22 


2.0 


2.4 








461 ± 92 






3.53 


3.0 


4.1 








245 + 49 






3.55 


3.2 


3.9 


239±29 


259±32* 










4.88 


4.5 


5.3 








128 + 26 






4.49 


4.0 


5.0 




144± 20 










5.73 


5.0 


6.5 




190±35 










7.87 


6.4 


9.5 


149+27 


151± 21 










8.3 


6.0 


10.9 












159±32 


12.0 


7.6 


15.4 










163+ 24 




23.7 


20.5 


28.5 


118±17 




107 + 11 








25.0 


16.5 


30.4 










108+ 16 




56.0 


38.6 


76.5 








700 ± 140 






60.0 


37.0 


82.8 










536 ± 80 




71.4 


55.1 


91.6 


1086 ± 256 




940 ± 188 








97.7 


68.6 


121.8 








3706 + 741 






100.0 


74.4 


130.6 










2928 ± 439 




155.9 


128.9 


184.0 


7315 ± 1632 




7900 + 1580 








147.9 


108.2 


181.7 








7545 ± 1509 






175.0 


140.0 


220.0 












7900 ± 800 


247.9 


174.7 


335.9 








6242 + 1248 







" The wavelengths reported in the first column are the nominal wavelengths of each instrument waveband. The minimum and maximum 
wavelengths reported in the second and third columns correspond to a system response (filter+detector) equal to 10% the maximum. In particular 
the system response functions were taken for Spitzer/IRAC from http://ssc.spitzer.caltech.edu/irac/spectraLresponse.html, for Spitzer/MIPS from 
http://ssc.spitzer.caltech.edu/mips/MIPSfiltsumm.txt, for COBE/DIRBE from http://lambda.gsfc.nasa.gov/data/cobe/dirbe/ancil/spec_resp/DIRBE 
_SYSTEM.SPECTRAL_RESPONSE_TABLE.ASC, for IRAS from http://irsa.ipac.caltech.edu/IRASdocs/exp.sup/ch2/tabC5.html, for 
MSX/SpiritHI from http://irsa.ipac.caltech.edu/data/MSX/docs/MSX_psc_es.pdf (Appendix A, Table A-3), for ISO/ISOPHOT from 
http://www.mpia.de/ISO/welcome.html. 

* Uncertainties in this table have been directly taken from the literature when available. For COBE measurements we assumed an error equal to 
20% of the measured flux in each band as suggested by Odenwald et al. (1998) and for IRAS measurements equal to 15% as suggested by Rice et 
al. (1988). For our own measurements the uncertainties have been derived as detailed in the text 

Table 3. Average values and standard deviation (in parenthesis) of lg(TIR/FUV), (FUV - ;'sdss) and A F uv in each annular region shown in Fig. 12 
(upper panel). The fifth and sixth columns show the linear regression parameters obtained fitting the data in Fig. 12 (bottom panel) with the relation 
lg(TIR/FUV) = a (FUV - ;'sdss) + b. The last column shows the number of regions used in each annulus. 



R(kpc) 


lg(TIR/FUV) 


(FUV - i SDSS ) 


^fuv 


a 


b 


N.Reg. 


1. 


2.13 (0.11) 


7.51 (0.07) 


1.20 (0.14) 






9 


3. 


2.15 (0.16) 


7.18(0.45) 


1.31 (0.16) 


0.38 (0.13) 


-0.55 (0.93) 


115 


5. 


2.12(0.28) 


6.68 (0.70) 


1.50 (0.30) 


0.41 (0.11) 


-0.62 (0.74) 


381 


7. 


2.05 (0.23) 


6.18(0.75) 


1.71 (0.40) 


0.32 (0.10) 


0.10(0.64) 


590 


9. 


2.03 (0.26) 


5.66 (0.71) 


2.06 (0.44) 


0.37 (0.15) 


-0.05 (0.88) 


879 


11. 


1.91 (0.24) 


4.96 (0.68) 


2.45 (0.47) 


0.42 (0.24) 


-0.19(1.21) 


1200 


13. 


1.81 (0.26) 


4.84 (0.69) 


2.37 (0.47) 


0.47 (0.27) 


-0.45 (1.35) 


1228 


15. 


1.61 (0.22) 


4.63 (0.64) 


2.17(0.28) 


0.35 (0.10) 


0.03 (0.47) 


1380 


17. 


1.59 (0.21) 


4.66 (0.80) 


2.06 (0.33) 


0.28 (0.09) 


0.29 (0.42) 


615 



